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Abstract 

We apply the effective field theoretic (EFT) approach to resum the large perturbative logarithms 
arising when partonic hard scattering cross sections are taken to the threshold limit. We consider 
deep inelastic scattering, Drell-Yan lepton pair production and the standard model Higgs produc- 
tion through gluon-gluon fusion via a heavy-top quark loop. We demonstrate the equivalence of the 
EFT approach with the more conventional, factorization-based methods to all logarithmic accura- 
cies and to all orders in perturbation theory. Specific EFT results are shown for the resummation 
up to next-to-next-to-next-to leading logarithmic accuracy for the above-mentioned processes. We 
emphasize the relative simplicity by which we derive most of the results and more importantly 
their clear physical origin. We find a new relation between the functions in the quark and 

gluon form factors and the matching coefficients in Drell-Yan and Higgs production, which may 
explain their universality believed to hold to all orders in perturbation theory. 
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I. INTRODUCTION 



Factorization theorems for inclusive hard scattering processes are our main tool by 
which we quantitatively analyze cross sections, with hadrons involved, when a generic hard 
scale is taken to infinity. Taking the Drell-Yan (DY) lepton-pair production as an ex- 
ample, it is well known that the cross section can be expressed as a convolution of a hard 
scattering coefficient function (or "Wilson coefficient" ) calculated perturbatively, and a non- 
perturbative, universal, parton distribution function (PDF) for each one of the incoming 
hadrons. Corrections to the factorized cross section scale as l/((5^)" where n > 1 up to 
some logarithmic ratios. However, it is also well known that fixed order, pQCD calculation 
of the Wilson coefficient yields singular distribution functions of the form 



where z = Q'^/s and s is the total momentum squared of the incoming partons. The 
"plus" distributions are defined in the usual way. The appearance of such distributions 
is a result of an emission of soft and/or coUinear gluons into the final state. When such 
distributions are Mellin transformed to the conjugate space, logarithms of the form In™" 
(m = 2k, 2k — 1, 0) show up where N = N expi^^E) is the conjugate variable of z and 7^ 
is the Euler constant. In the limit 2; — 1 or, equivalently, large N, fixed order perturbative 
calculation cannot be reliably trusted and an all order resummation of the large logarithms 
is needed. This is what is generically meant by "threshold resummation" . This notion has 
evolved during the last 20 years into one of the most studied and highly developed subjects 
within perurbative quantum chromodynamics (pQCD). Earlier studies 0,13] supplied a sound 
and rigorous (although complicated) treatment to perform such resummation. In both of 
these works, resummation is performed after establishing some sort of factorized cross section 
into well-defined quantities (at the operator level) that capture the physics at the hard, jet 
and soft scales. An integral transformation to the conjugate space is then applied in order 
to de-convolute the various terms in the cross section. Then, in the conjugate space, energy 
evolution equations are solved and the exponentials thus obtained contain the resummed 
large logarithms. Thus the perturbative expansion is put under control and the contributions 
obtained from the yet uncalculated higher orders in reduce the theoretical uncertainty 
inherent in any fixed order calculation. Thereby better phenomenological studies can be 
carried out and better agreement with experimental data is usually obtained. More recent 
studies have further developed and refined this topic 0, IE S S H 01 • 

In this paper we adopt the effective field theory approach (EFT) to resum the threshold 
large logarithms. This approach was first applied to deep inelastic (non-singlet) structure 
function in the limit x ^ 1 where x is the par tonic Bjorken variable jl^]. Later on it was 



applied to the DY process in the limit z ^ 1 In both cases, resummation was performed 



up to next-to-leading logarithms (NLL). The implementation of the EFT methodology to re- 
sum threshold logarithms is made more concrete due to the recently developed "soft collinear 
effective theory" (SCET) 

The SCET describes interactions between soft and collinear partons. It is the most 
appropriate framework to calculate contributions from the soft-coUinear limit of the full 
QCD calculations (which is more commonly known as the "soft limit"). Therefore any 
perturbative calculations within SCET has to reproduce the same results as the full QCD in 
that limit. To 0{as), this has been verified explicitly for DIS and DY [lol llll|. This is also 
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the case when one considers distributions at small transverse momentum 3 3 • Moreover, 
it was also shown that the one loop diagrams (the form factor type of diagrams) calculated 
in SCET have the same infrared (IR) pole structure as the full QCD calculation 0, 0|. 
[The result in for the collinear diagrams involve mixed poles of IR and ultraviolet 
(UV) divergences. This would be treated by applying the "zero-bin" subtraction jl7l|]. 
These observations have to be valid to higher orders in the strong coupling. This allows 
us to extract the relevant quantities needed to perform resummation from the full QCD 
calculations as we shall see below. 

The EFT resummation program described here is conceptually simple and has been 
explained in detail in The starting point (and again considering DY as an example) is 
the collinearly factorized inclusive cross section in moment space [19| 

aN = ao- Gn{Q) ■ q{Q, N) ■ q{Q, N), (2) 

where cxo is the Born level cross section, q{Q, N) is the PDF of partons in hadrons and, 

Gn{Q) = |C(«,(Q2))|'e^i(«/'^^'"^(«')) X 7Wiv(«.(Q'))e'^(^/^^'"^(^?))e'^(«/^^'"=(«». (3) 

Explicit expressions for the various contributions in G^v will be given below; however, we 
want to comment on their physical origin. C(q;s((5^)) contains the non-logarithmic contri- 
bution of the purely virtual diagrams and the first exponent Ji contains all the logarithms 
originating from the same type of diagrams. Both quantities are obtained from the matching 
procedure at the scale Q, and the running between Q and the intermediate scale /i/. This, 
in turn, is controlled by the anomalous dimension of the EFT current to be denoted by 71. 

The intermediate scale /// shows up when real gluons are emitted so one must consider 
the cross section with real gluon emissions. The result will have both soft and collinear 
divergences. When taking into account the IR poles from the virtual diagrams (in the EFT 
approach this is done by taking into account the contribution from the counterterms of 
the effective operators), the total contribution will contain only collinear divergences to be 
absorbed into a product of two PDFs. The conclusion is that the matching procedure at 
the intermediate scale is guaranteed to work to all orders in perturbation theory, following 
the factorization theorem, as long as the EFT used generates the full QCD results in the 
appropriate kinematical limit and one gets the matching coefficient Ai]^{as{fij)) which by 
definition is finite in the non-regulated theory. This quantity has to be free of any logarithms. 
I2 collects all the the logarithms that are due to the evolution of the PDF between fij and 
the factorization scale fip- This is controlled by the anomalous dimension of the PDF to be 
denoted by 72 . I3 encodes all the contributions due to the running of the coupling constant 
between the matching scales {Q and ni) and the final factorization scale fip- All the large 
logarithms appear only in the exponents and the term |C(as(<5^))P-^Jv(«s(Q^)) is free of 
any large logarithms. In Eq. Q we have chosen fip = Q for simplicity. 

The above formalism must be contrasted with the more conventional, factorization based 
one. It will be shown that the EFT approach is equivalent to the other approaches and to all 
logarithmic accuracies. We will derive all the known ingredients needed to perform threshold 
resummation to next-to-next-to-next-to leading logarithmic accuracies (N^LL) for DIS non- 
singlet structure function, DY process and the closely related Standard Model (SM) Higgs 
production through gluon-gluon fusion into a top quark loop. Moreover, the integrations 
in Eq. (jS)) are very easy to perform and we have carried out the integrations up to g^^^ 
that resums the NNLL. This calculation is to be compared with the ones explained in, e.g.. 
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Appendices A of Refs. 2^- In this paper we use dimensional regularization in d = 4 — 2£ 
to regulate both the UV and the IR divergences and we utilize the MS scheme throughout. 

This paper is organized as follows. In Sec. II we derive the anomalous dimension of the 
quark and gluon effective currents up to and write down the matching coefficients 

at the scale up to O^a"^). In Sec. Ill we obtain the matching coefficients at fij to C(o;^) 
and give our final expression for the resummed coefficient function G^. There we also 
comment on the universality of the functions that enter the quark and gluon anomalous 
dimensions of the effective operators. In Sec. IV we compare the EFT approach with the 
conventional one and derive our main result that establishes the full equivalence of the 
two approaches. From that relation we obtain the recently calculated D(^q ioi DY and 

(3) 

Higgs production and Bq for DIS. We carry out the integration in the resummed coefficient 
function to illustrate the simplicity of the EFT results and obtain the well-known functions 
g^^\^){q,g) for i = 1,2, 3. Our conclusions are presented in Sec.V. In the Appendix we write 
down explicit expressions for soft and virtual limit in full QCD for all the processes we 
consider up to O^a^) in z space and in the conjugate space (for large moments). 



II. ANOMALOUS DIMENSION AND MATCHING COEFFICIENTS FOR EF- 
FECTIVE CURRENTS 

The EFT approach for resummation starts from calculating the contributions at scale 
Q^. Technically this is done by matching the full QCD theory currents to the EFT currents 
at the scale by considering the purely virtual diagrams in the full theory. By doing this, 
we integrate out the hard modes of virtualities of order Q^. The matching of the currents 
can be expressed as an operator expansion 

Jqcd = CiQ^fi^, Jeff(/i) + (4) 

where C is the matching coefficient, /i is the factorization or renormalization scale of the 
effective current and ellipses denote higher-dimensional currents which will be ignored in 
this work. We will consider the quark vector current J'^ = 'ip^^^ij) for DIS and DY cases and 
the gluon scalar current J = G'^'^G^y for Higgs production in hadron colliders. 

The anomalous dimensions of the effective currents that control the running (with ix) are 
defined as 

din Jeff /ex 

7i(/i) = (5) 

If the matrix elements of the currents in full QCD are independent of the factorization scale, 
such as quark vector and axial vector currents, the same anomalous dimensions are obtained 
from the matching coefficients of the effective currents 

, , dlnC 

7i(/^) = /i^. (6) 

The anomalous dimension is a function of both Q^Z/i^ and Usi^"^)- In fact, it can be shown 
that it is a linear function of lnQ^//i^ to all orders in perturbation theory [lo| : 

7i = A(a,)lnQV/^' + 5i(«s), (7) 

where A and Bi have expansions in = Q;5(/i^)/47r: A = '^s^*'*^ ^^'^ -^i — Xli '^s-^i*^ 
and as is the renormalized coupling constant. 
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To obtain the anomalous dimensions and the matching coefficients, we consider the sim- 
plest matrix element of the full QCD currents between on-shell massless quark and gluon 
states. They are just the on-shell form factors F. Since they are "physical" observables, 
there are no UV divergences, but there are IR ones. To all orders in a^, we can write 



(8) 



where S contains only infrared poles in dimensional regularization (i.e., no finite terms). S 
can be regarded as the matrix element of the effective current Jes after renormalization has 
already been performed. In the effective theory, Feynman diagrams for S have vanishing 
contributions in dimensional regularization because there are no scales in the integrals. This 
can be regarded as the result of cancellation of IR and UV poles. As such, the IR poles in 
S may be treated as UV poles for the purpose of calculating the anomalous dimension 

din 5* . . 

71 (/i) = -/^— (9) 



Since C does not contain any pole part, we can also write 

d\nF 



7i(/^) 



(10) 



pole part 



Therefore, the perturbative results for F up to any loop order can be used to calculate the 
anomalous dimensions to the same order. 

The best way to see the physical content of the form factor is to consider a resummed 
form m H 13 



In F{as 



^ ( Kiasif^), e) + G{1, asi^fx, e), e) + 



yv4(as(A/i,e; 



where A is the anomalous dimension of the K and G functions. 



dG 



,dK 



/i 



(12) 



and is in fact the same A as in Eq. ((Tj). K contains only the IR poles, and therefore, the 
whole ii'-function can be constructed from the perturbative expansion A = Yli'^l^^^^- The 
function G contains only the hard contribution, and has a perturbative expansion 



G(l,a„e) = 5^a:G»(e) 



(13) 



Thus In F can be expressed entirely in terms of G'-*-' and A*^*) . 

The anomalous dimensions Ag for the quark vector current (DIS and DY) and Ag for the 
gluon scalar current (Higgs production) have been calculated up to 0{al) |24| . 



A 



A 



A 



(1) 

(2) 
(3) 

{1:9) 



16C, 



+CaNf 



245 



In. 



24 
209 



67 11 

y C2 + y Cs 



10, 7, 

\ C2 Cs 

108 9 ^ 3^ 



11.2 
5C2 



+ CfNf{-- + 2(, 



^^'-27 



(14) 
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where C(q^g) = Cp for the quark and Ca for the gluon. In this sense A is universal. 

The expansion coefficients for the G function have been obtained up to 3-loops from 
exphcit calculations of the quark and gluon form factors 25j : 

^Z) = - 35,/32) + C) + + /5o - f3o^\'jJ , (15) 



'{1,9) ^"^(m) 



where 6g is zero for quark and 1 for gluon. The B2's are the coefficients in front of the delta 
function 5(a; — 1) in the Altarelli-Parisi splitting function and have been calculated to the 
third order l2j 



Bfl = ACfCa + yC2 - 3C3) - ^CpNp + + 4C| - 3C2 + 6(3 



Bi'l = IQCaCpNf ( - - — C2 + + -Cs 

2,q A t t S2 20^2 ^gS3 

^^(^r (^^^ ^ 205 247 211 15 
+IQCaC, ^— + C2C3 - ^C2 - -^C2 + ^C3 + yC5 

-IQC^C, - — C2 + gC2 + yC3 - 



^ ^V144 27^ 9^^ 
-lec^iv^ - -C2 - ^C2 + yC3 

+16CI g - 2C2C3 + ^C2 + f CI + y C3 - 15C5 ) , (16) 

for quarks, and 

B^l = -Ca - -iV^, 

= 4CaNp (-1) + 4Cl f § + 3(3] + ^CfNf ^ ^ 



3 / V3 7 V 2 

288 J ^288 ^ V288 6 12 



= WCaCfNf ( -^^^ ) + IGCaN^f^ - IGClNp ( + ^(2 + ^Cl + ^Cs 



+16Ci - C2C3 + ^C2 + ^Cl + f C3 - 5C. 
+16C^iV|^ + WC^Nf^, (17) 
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for gluons. The universal functions f{q,g) are given by 



(1) 
(9,9) 



f{2) 
J I 



(q,9) 



A3) 



0, 

C{q,g)CA 
^{1,9)^ A 



808 22 
136781 12650 



729 

'^^(q,9)^aNf 



C2 - 28C3 
C2 



112 8 



1361 352 . 176,, 

-^Cs + — Cl + — C2C3 + 192C5 



+4C2 + 



304 



9 



C3 + 



11842 2828 , 

\ C2 

729 81 ^ 

32 

5 



728 



96 



27 5 



+ ^(q,9)^P^P 



1771 



2080 40, 112, 

C2 H Ca 

729 27 27 



(18) 



The tilted functions in Eq. ()15|1 are not given here since they do not contribute to the 
anomalous dimension as their contribution to the form factors are canceled (they can be 
found in 111). 

Finally, the anomalous dimension of the effective currents can be expressed in terms of 
the A and G functions. If one writes 71 = al'yi \ then 



where 



B 



1,(9,9) 



-25, 



U<1,9) 



(9,5)- 



and the QCD /5-function is given by 



I3{as 



dhias 

din yU^ 



(3oas + I3ial + 



(19) 
(20) 

(21) 



with (3o = 11Ca/3 — 2Nf/3. The above expression for 71 might work to all orders in per- 
turbation theory. In the gluon case, the last term is present when the anomalous dimension 
is defined in terms of the matching coefficient Cg and is absent when it is defined in term of 
the effective current. 

The anomalous dimensions could also be calculated from the matching coefficient 
C(Q^//U^, as(/i^)) extracted from known results of the form factors. First, we take the 
logarithm of Eq. (jH)) , 

InF = InCiQ^fi^) + lnS(gV/i^ 1/e). (22) 

Then we separate out the poles from the form factor logarithms, which belong to the 
S'(Q^//i^, 1/e). The finite part left over is just the logarithm of the matching coefficient 
In C to any desired order. So, eventually, we will get the following result for the anomalous 
dimension valid to arbitrary order in a^. 



7i 



d 



{\nF\ 



finite part} ; 



(23) 



where the form factor in the above equation has been renormalized (including coupling 
constant renormalization). Using the above equation, we have calculated the anomalous di- 
mension for the quark and gluon currents in the effective theory up to three-loop order, and 
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they are exactly the same as Eq. (fT^. We should point out here that the anomalous dimen- 
sion of the quark current is the same for the scattering case (DIS) and for the annihilation 
case (DY). 

To calculate the matching coefficients, C((5^//i^, ^^(/i^)) = J2i (^lif^'^)'^^^KQ^ / ^'^) fo'^ DIS, 
DY and Higgs production, we need the expressions for the quar k form factor (space-like case 



and time-like case) 129] , and for SM Higgs production j27|, |28|, |29| up to the same order. It 
should be noted that for DY and Higgs cases, C^*^ contains imaginary parts that need be 
taken into account. For our purposes, it is enough to keep the imaginary part for C^^-* only. 
Normalizing C^^^ to 1 we find for DIS, 



cl 



-In^ 



+ 3 In 



-8+C2 



ln^p_ -31nP^ +8-C: 



+ ( - - I2C2 + 24C3 



+CfNf 
4085 



2, 3 
— In"* 
9 



+ 



23^ 

\ C2 

324 9 ^ 



+CfCa 
f 254:5 



11 



In^ 



+ l^ln^ 
9 



1 44 
- + 29C2 - 3OC3 - -('2 



+ 2C2 - 



Q: 



233\ 

Is) 



22 



26C: 



In^ 

51157 



209 4 
^ + 3^^ 



In 



(24) 



337, 313 

C2 H 

18 ^ 9 



^ 44^2 



648 18 ' 9 ' 5 
The logarithms in the above result have been presented in 3|. For DY we can simply get 



the Cq^ by replacing each In (^^j in C'dis "with In (^^j ~ "^^r. This is just a result of the 
fact that the time-like quark form factor can be obtained from the space-like one by analytic 
continuation. For the Higgs production we set Mjj = and we get 



^2 



+ 



2 
80 
27 



+ 7C2 + 2mMn( ^ 



^C2 



2C3 In 



67 

I7C2 

9 



+ 



5105 335 

+ 



143 



125 



+CaNf 
916 



9 



H In^ 

9 



162 



6 



C2 - ^C3 + —rC 



9 



52 16, 

\ C2 

27 3 ^ 



10 



81 

+CfNf 



25 46/ 
y C2 - yCa 



21n^ 



/i2 



67 



+ 8C3 



(25) 



The logarithms in C*^*-* will be needed later on to show that the matching coefficients at fij 
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are free of any logarithms, and we have not included the imaginary part of Cg since it does 
not contribute to the accuracy in which we are interested. 

Using Eq. (0) we can write down the solution of the renormalization group equations for 
DY and Higgs, 



h{Q,IJ'i) 



(26) 



C{q,g){Q^/^^^i,as{^^j)) = C(q,g)(l,a4Q^)) exp 
where 

^1 = - / — 

7i,(g,9) = 7i,(<?,9) - 2i5gA-i- (27) 

The C(l,as((5^)) = C(as((5^)) is just the non-logarithmic part of C{Q'^/fi^,as{^^)). For 
the Higgs case the last relation is a result of the fi dependence of C^{fi) which enters into 
the effective lagrangian that one obtains after integrating out the top quark (see, e.g., (27^). 
This /i-dependence is governed by anomalous dimension which we denote by jt following 
the notation of There it was shown that 

7T = a.[-2/5o]+a'[-4/3i] , (28) 

so the conclusion is that the only effect of this anomalous dimension, when combined with 
anomalous dimension of the matching coefficient at the scale is to cancel the terms in 
7i for the Higgs case. For DIS we replace Cg with Cdis which runs with the same 7i,g. 

In Eq. ()26p we encounter the first of three exponentials. The other two will be obtained 
below. Since fij will later be identified with Q"^ /'W where p = 1 for DIS and p = 2 for DY and 
Higgs cases, it is clear that the exponential includes large logarithms of the form mentioned 
in the introduction. We again stress the fact that C(as((5^)) (for all three processes) and 
'Ji,{q,g) are completely determined to a given O^a^) by the knowledge of the form-factor 
calculation up to the same order. 



III. MATCHING COEFFICIENTS AT fii AND THE RESUMMED COEFFICIENT 
FUNCTIONS 

In this section we show how to extract the matching coefficients at the intermediate scale 
to 0{al) for DIS, DY and Higgs production from the known calculations of full QCD, and 
to obtain resummed expressions for the coefficient functions. Since we are interested in the 
threshold region, we need to consider only the partonic channels that give rise to the singular 
contributions in the limit z ^ 1, i.e, 5(1 — z) and the "plus" distributions, 'Di{z), where 



VAz) = 



In'fl - z) 



(29) 



For DIS, DY and the Higgs processes, these channels are: q + 'J* q, q + q 'j* and 
g + g ^ H, respectively. 

To the accuracy we are interested in, the 0(Q;f) cross section from soft contributions are 
needed. The full QCD calculations for cross sections can be found in Refs. [26|] for DY, in 
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Refs. j28|, |29| for the Higgs production and in Refs. [30|, |31| for DIS. The result in the soft 
hmit can be written as 

G^'+'^z) = J2<if^^)G^'^'^"^^\^) (30) 

i 

which contains both soft and virtual contributions, and where G{z) is the inverse Mellin 
transform of Gat in Eq. (j2)). Explicit expressions for G^'^^+^H^ with i = 1,2 can be found 
in Ref. [s^ for DIS, in Ref. [s^ for DY and in Ref. [s^ 3^ for the Higgs production. 



[G^'^\z) = 5{1 — z)]. Using the following well-known Mellin transforms of 'Di(z) in the large 
N limit 

Vo(N) = -IniV, 

V^(N) = -lln^AT-CalniV-^Cs, 

V,{N) = lln^iV+^C2ln'iV + 2C3lniV+^C2, (31) 

we get the G^^^'^^~^^\N), i = 0, 1,2. Explicit expressions are given in the Appendix. As we 
have already mentioned, the SCET is supposed to reproduce the same results. 

To get the matching coefficient at the intermediate scale /i/, M.^ = X]j '^s-^!v''> need 
to factorize the virtual contribution from the following relation. 



(s+v) [Qjrf 



-|2 



2 /n2 



G'^' ( ^, N, a.ifi'yj = |C «s(/^') j f X -^^ ^' "«(/^') ) • (32) 



The content of this formula is simple: The finite part of the partonic cross section G^^^'' 
comes from both the purely virtual, form-factor type of Feynman diagrams, which are in- 
cluded in |Gp, and from diagrams with at least one real gluon emitted into the final state, 
which are included in Ai^. However, there is a different way to look at it. The right-hand 
side is just the result of a two-step matching of the product of two full QCD currents where 
at each step we collect the relevant contribution to the cross section. The first step accounts 
for |Gp and the second one gives rise to Ai^. It should be noted that multiple matching 
procedure, as the one performed here, results in a multiplicative matching coefficients. We 
also mention that the above equation could formally be proved, inductively in a^, by con- 
sidering the cross section within the effective theory itself and relating it to the full QCD 
calculation in the soft limit. 

Expanding the above equation to the third order, one gets 



^{2),s+v ^ |^(i)|2 ^ 2Re[G(2)] + 2Re[G(i)]^£) + (33) 
gS?)'^+^ = 2Re[G(i)G(2)*] + 2Re[G(=')] + |G«|2A<ii) 
+2Re[G(^)]^g^ + 2Re[G(2)];wi;^ + 

The above factorization is consistent with that considered in j^^. We get the following 
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result for DIS, 



M 



(1) 

Af,DIS 



•'^liV,DIS 



Cf [2L^ + 3L + 7 - 4C2] 
37 



cl 



+CfCa 



2\} + 61^ + 
22. 



9 



2 

367 

Is" 



- 8C2 J L2 + 

- 4C2I - 



45 



- 24C2 + 24C3 L 



3155 22^ ^ , ^ 
.__ + _C2 + 40C3)L 



9 9 



205 97 

— -yC2-6C3 



4, 247 

-C2 

3 27 



122 



L 

CfCa 



53129 155 



37 



648 



-7rC2 - 18C3 - ^C2 

6 5 



4057 13 , 

^ C2 

324 3 ^ 



(34) 



where L = In The above resuh has also been obtained in where an explicit two-loop 



calculation of a suitably defined jet function was performed j43[. For DY, we get 

Cf + 2C2] , 



[2L^ + 2C2Y + CaCf 



22 



134 



4C2 L2 + 



808 



28C3 L 



+CfCa 
+CfNf 



9 



20, 
"9 



2428 67 



112, 

27 
22 



yC2 - y C3 - 12(2^ 



328 10 , 4 , 

C2 + -Cs 

81 9 ^ 9^ 



(35) 



where L = In ■ And finally, for the Higgs case, we have 



= Ca[2L' + 2C2\ 



di -)[2L' + 2(2? + CaCa 



2 

-CaNf 
+CaCa 
+CaNf 



22^3 
9 



134 



4C2 L 



808 



9 9 



2428 . 67 
9 



27 
22 



+ ^C2 - ^CS - I2C2 



328 10, 4, 

C2 + -Cs 

81 9 9 



For all three processes, we have G^^^ = 1. 

From the above results it is clear that the logarithms L and L vanish when we set: 
/i^ = /if = Q"^ I'W . Of course, this has to be the case as the matching coefficients should be 



(36) 
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logarithmically free, and we can write 



Mn (j^, iV, = Mn (\n , «.(/^')) , (37) 



and for fi"^ = fij = ^ we have 



( In ( =^ ) , asifx') ] = MNiasifx'j)). (38) 



Q 

These observations are valid to all orders in perturbation theoryfl^ and they lead to a 
strong constraint on the anomalous dimensions of the effective operators on both sides of 
the matching scale. Another interesting feature emerges from the results of the DY and Higgs 
cases, M'-^lg and A^g^, i = l,2: One can simply get the latter from the former by replacing 
the overall factor Cp with Ca in the non-Abelian part. The Abelian part exponentiates and 
hence all occurrence of Cp shall be replaced by Ca- In this sense, the matching coefficients 
seem to be universal. This could be argued based on that in the soft gluon hmit, only the 
color charges of annihilating quarks and gluons are relevant. 

Following the same steps as we did after the first stage matching at Q^, we need now 
to consider the running of the effective operators that were used to perform the matching 
at fij. However at and below the scale /// they are just the conventional PDFs taken to 
the limit z ^ 1. As such, the running of the effective operators (the PDFs) is governed 
by the well-known DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parizi) evolution equation 
with anomalous dimension 

72^fe,) = %.)lniv'-252,fe,), (39) 
where and -82,(^,9) are given in Eqs. (fT^ and (fT^ . We include the running effects in 

h = 2 ^72,(.,), (40) 

where hf is the factorization scale for parton distributions. 

The resummed factorization coefficient functions for DY and Higgs are 

GN,i,,dQ) = iq,,,)(«s(Q))re^^(«'^^) X MN,i,,,)Mf^i)y'^''"''\ (41) 

where we have omitted C| for Higgs production. [The definition of Ji and I2 differs by a 
minus sign from Ref. Anticipating the discussion of the next section, we will set the 

factorization scale hf = Q- The above equation can be brought into an equivalent form by 
exploiting the running of from to Q in A4N,{q,g){c(s{fJ'i))] 

MN,ig,9)i(^s{f^j)) = MN,ig,9)MQ^)) ewlh], (42) 

where 

h = -2 r^Ai?(,,,), (43) 
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where 



The last two equations are also true for the DIS case. Thus we write 

G;v(Q)=-^(«.(Q))e^(''"^(«)\ (45) 

where T = \C(^g^g){as{Q))\'^A4{g^g){as{Q)) depends only on as{Q). The subscript of has 
been omitted since there is not any large logarithmic dependence in the matching coefficients. 
/ = /i + /2 + /s is a function of X = Po In Nas{Q) and as{Q) with all leading and sub- leading 
large logarithms resummed. 

Since the cross section ctat in Eq. (2) is independent of the intermediate scale /i/, then 
from Eq. ()41|) and the definitions of 71 and 72 we get the following relation for DY and 
Higgs; 



d\llMN,{q,g){(^s{fJ''^),l>) 

from which we get 



[272 - 27i](,,,) = 2[AL + /](,,,), (46) 



din fx 



= 2/(,,,)(a.(/z?)), = ! (47) 



where are given in Eq. (fT^ and /(g^^) are given in Eq. (fTHj). The last equation sheds 

light on the physical meaning of the functions f{q,g)'- It is the anomalous dimension of 
the matching coefficient A4 evaluated at the intermediate scale fij. Here we see that the 
universality of these functions could be explained by the fact that ^A(q,g) are themselves 
universal. 

The last equation also shows the same appears in the logarithmic parts of 7i,(g,g) 

and 72,(g,g), because otherwise the logarithms at /i/ do not cancel in A4n- 

For DIS a similar analysis is performed, however, we have to consider only one-half of I2 
in Eq. (|4ip since we match onto a single PDF. With this we get 

from which we obtain at the intermediate scale, 



dliaMN,Bis{ois{fi),L) 



dlnn 

Here there is an extra contribution from 



= 2[i?2 + /],(a,(/x?)), fij = -^. (49) 



IV. COMPARISON WITH THE TRADITIONAL APPROACH AND EXPLICIT 
RESULTS TO N^LL ORDER 

In this section we will illustrate the equivalence of the EFT approach and the traditional 
one which relies on the refactorization of hard processes as we mentioned in the introduction. 
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The renormalon problem in the later approach arises from doing resummation uniformly 
for all moments, which will necessarily encounter small scale Q/N^ at fixed Q when is 
sufficiently large. The EFT approach avoids that by short-cutting the steps when this scale 
becomes of order Aqcd- We will start by showing this first for DY and Higgs production, 
then will turn to the DIS case. In the last subsection, we give the explicit form of the 
relevant integrals obtained in the EFT approach. 



A. Drell-Yan and Higgs 

One of the well-known forms used to express the coefficient function for DY and Higgs 
in moment space is the following js^]: 



GNiQ') = goMQ'))e'^AC{as{Q')), (50) 

where we have normalized the Born term to 1. The Qq has a conventional expansion form: 
go = J2i O'^Qoi- [III this subsection, we omitted the subscript q and g, intended for DY 
and Higgs production.] The term AC has the only role of cancelling the non-logarithmic 
contributions that appear in the exponent. These contributions arise from the various (- 
terms in the Mellin transform of the "plus" distributions. The Sudakov exponential term 
I A is given by 



'A 



dz- 



1 - z 



-±-Aia^fi')) + Diasiil-zyQ')) 

|2 yU 



(51) 



where, as already mentioned, we set /i^^ = Q^. As noted above, Ja contains both a logarith- 
mic and non-logarithmic contribution. The quantities, go, A and D have the usual expansion 
in as and they are already known up to 0(af) |20|. The A is identical to the logarithmic 
coefficient in 71 and 72. It is our aim to relate these quantities with those that appear in 
Gn of Eq. (jH}. For this we follow the procedure outlined in Appendices A, B and C of |l9l |. 
The integral in can be rewritten in terms of the already defined /i, I2 and J3, 

I = h + h + h = lA + In AC(a,(g2)), (52) 

where the coefficient function AC does not depend on fij ^ Q/N. 
To prove the above relation, we first use following expansion; 

where the F function is related to the usual gamma function, 

r(.-J^]..-rJJ^](J^)\ (54) 



dlnNj \dlnNj V^lnA^, 

where the first parenthesis in the right-hand side is the argument of the F2 function, and 

r2(6) = ^[1 - e-^-F(l - e)] = -1(2 - ^Cae - ^^e' + 0{e% (55) 
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In Eq. (jH!^ we used (d/d\nN)f(\nN) = {d / d In N) f (In N) for an arbitrary function /. 
After some algebra, /a can be expressed as 



'A 



-r 1 



d 



AK(/.2))ln^ + iD(/.2) 



(56) 



The double derivative from F acting on the curly bracket above gives a contribution 



d 



d\nN 



d 



dluN 



D(a,(QViv'))-4A(a,(QViv')) 



(57) 



To compare /a with the exponent / = /i + /2 + /s, we express the latter in the form 



h+l2 + h 



A{as{ii^)) In ^ + (fii + AB + 2^2) 



A(a,(/i2))lnAr 



(58) 



Matching the two integrals, we get 



J- 2 



(91nA^ 



d 



d\nN 



=D{a,{Q'/N'))-AA{a,{Q^/N')) 



.^1 ,,2 



D(«,(/i^))+lnAC(«,(g^)). 



The above equation can be solved by perturbative expansion in «<,. 

If the equality given in Eq. holds to all values of A^, then for = 1 we get 



(59) 



lnAC(«,(g2)) = -r2(9„J 9„,Z}(a,(gVArO) -4A(a,(gVAr')) 



(60) 



where we follow and replace the derivative d/d\\iN with das where 



9„ 



d\n.ij? das 



_d_ 



and, hence, 



,91nA^ 

where / is an arbitrary function. 



% ] /(a,(gVA^')) = daJ{aM/N\ 



(61) 



(62) 
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Applying one more d/dlnN = da^ on both sides of Eq. we get our master relation 

2(5i + AB + 252)(a.(/i')) = Z^(«.(^')) + 9,j2(9„J [AA - d^M (63) 

which can easily be solved for D^^^ order by order in a^. As an example, let us expand both 
sides up to O^a^). First, we work out the expansion of the AB term. From Eq. PHj). we 
get 

il,9) 



ABf^ , 
ABf^ , 



ABf^ , 

{1,9) 



^<9) = 0' 






2 M^'^^ — 










-Pi 







(2) 



Noticing that B^^^^^ + 2S«^^^) 



A^,(<?,9)- 

{1,9) 



f^j' ) and using the expansion of we get D^*-' 



(64) 



(65) 



D 



(0) 
^(2) 



D 



D 



(3) 
(9,9) 

(4) 
(9,9) 



D 



(1) 



0, 



(9,9) 

-<9)+2AC)+W2A«), 



-2/< 
8C 
"3 



.(3) 
'(9,9) 

(4) 



(1) , /i(2) , 32^ .2.(1) 



(4) 



(9,9) ^ "'-2f^U^(q,g) -r g sJMO--(g,g)) 
(3) , QA/? /I (2) , qoA,q2 .(2) 



'(9,9) - —(9,9) - 12C2/5oA-^) + 8C2/3iA- ) + 32C3/3o^A;^^^) 



(66) 



Thus, apart from the coupling-constant running effects, D is essentially —2/ = 2Bi + 4i?2. 

From the last two equations we see that in order to get D^''\ the only same order infor- 
mation needed is f^^\ All the quantities needed to calculate D^"^^ and D^^^ are known and 
we get 



^(2) 
(9,9) 



C'(9,9) i 



101 11 7 \ /14 2^ 

^ C2 + -Cs + A^F C2 

27 3 ^ 2^V V 27 3^ 



(67) 



(3) 
(9,9) 



C(g,9)Cl 



594058 98224 , 40144 , 

+ — + ^^CS 



2992 



352 



+ C(g,g)CANF 



729 
"125252 
729 
3422 



11 
29392 
81 



C2 



27 

2480 
" ~9~ 



. , -CI - ^(2(3 - 384C5 

15 6 



. 736^2 
C3 + — C2 



15 



27 ^ 9^5 



+^(9,9)^1- 



3712 640 320 

1 C2 H G 

729 27 ^ 27^^ 



(68) 



16 



where C(g g) = Cp for the DY case and Ca for the Higgs case. The above results agree with 
the recent calculation in 41 1. The result for the Higgs production has already been 

reported on in [Ts^ . 



The non-logarithmic contribution J^{q^g){Q'^) = Xli g) = \C{Q'^)\'^Mn{Q'^) can be 



calculated from the already-known results for C'^^JgAQ'^) and J^^^', Jas{Q'^)), or we can 



simply read them from the well-known results for through Eq. ()32|) and Eq. ()33p : 



16Cf{C2 - 1) 
"511 



4 

+CfCa 
+CfNf 



198C2 
1535 



6OC3 



552 



376 
12" + —^^ 



-a 



604 92 2' 



127 
~6~ 



64 8/ 
y C2 + gC3 



(69) 



for DY lepton-pair production. For the Higgs case, we have 

I6C2CA, 



17(1) 
■^9 



(2) 



^2 



1535 376, 

\ C2 

12 3 ^ 



604 



9 



92 



80 160 88/ 



^CfNf 



67 



^ + I6C3 



(70) 



The above results agree with the goi and go2 in j38[. The ^f terms in the results of [38| are 
due to the use of iV instead of as in our case. It is very simple to also reproduce these 
terms. We also notice that their results for the g^i do not include the contributions from the 
non- logarithmic terms in J^. 



B. DIS 



For the DIS case there are essentially two major differences. The first is that the D 
term in is zero to all orders in jHliO]- The second one comes from the "jet function" 
which encodes the effects of coUinear gluon emission from the outgoing parton. So for DIS, 
the traditional approach yields the following expression for the exponent in the coefficient 
function Gn{Q'^), 



Ims = I dz — 



z 



(1-^)Q2 



dij/ 



A,{a,{^i^)) + B,{as{{l-z)Q^)) 



(71) 



where again we set 
introduced earlier. 



We have used B here so that it will not be confused with B^s 
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We now follow the same procedure as for the DY case, rewriting 



Idis = -r 1 



+ 



d 



I 

JQ 



d\nNj [Jq2/n /^^ 
2 A,(a,(/i2))ln]V}>, 



A,(a,(/x^))ln^+B,(/x^) 
A* 



'Q2 

On the other hand, our result for DIS reads 



7l + 72 + ^3 



o — //2 



A,(q;,(//2)) In ^ + (Si,, + ASdis + ^2,,) 



+ 



A,(a,(/i2))lniV 



Matching the two results above, and noting that 

d 



1 



we get the final relation between EFT and traditional approaches for the DIS case; 



from which we can solve for B^q . Up to third order we have 
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Therefore, apart from running effects, Bq is essentially — -82,9 — fq- More explicitly, we get 



-_ 

q 

BP = Cl 



-- + I2C2-24C3 



CpCa 



3155 44, 

+ — C2 + 4OC3 

54 3 S3 



247 8 



27 



;C2 



Bf = Cl 



29 

y 



288 



---- I8C2 - 68C3 - ~-^c,l + 32C2C3 + 240C5 



+C\Cf 
+CINf 



712 272 

-46 + 287C2 - ^Ca - ^Cl - I6C2C3 - I2OC5 
3 5 



599375 32126, 21032, 

+ ^;^C2 + ^^C3 



729 
5501 



32 

5OC2 + — C3 

54 ^ 9 ^ 



+CaCfNf 



160906 9920 



729 



-C2 



27 

+ CfNI 
776 



652 2 176 
^C2 - — 



C2C3 - 232C5 



8714 232 , 32 , 

\ C2 Cs 

729 27 ^ 27^ 



208 



-^C3 + ^C2 



(77) 



Those results agree with the ones in Ref. [20]. Similar to the case of DY and Higgs, we get 
after simple calculation 



(1) 



DIS 

(2) 



DIS 



16Cf(-9 - 2C2 
331 111 
2 



^ + ^(2 - 66C3 + 7C2 



4 
-( 

5 



+CfCa 
+CfNf 



5465 1139 



72 



464 



C2 + ^C3 + ^ci 



51 



Again these results agree with and 



457 85 , 4 , 

\ C2 + -Cs 

36 9 ^ 9^ 

,DIS 



(78) 



C. Drell-Yan Coefficient Function Using DIS Parton Distributions 



If one calculates the Drell-Yan coefficient function in terms of the DIS parton distribu- 
tions, one has 



GN,q/ G%r)is 



dz- 



1 - z 



(1-2)Q2 



+Dq{as{{l - zfQ^) - 2Bq{{l - z)Q^)\ 



(79) 



We have seen from the last two subsections that if one ignores the running effects, Dq ~ 
2Bi + 4i?2 and Bq ^ Bi + B2. Hence the last two terms in the above equation is just ~ 2i?2 
in EFT, negative of the coefficient in front of 5(1 — x) in the DGLAP splitting function. 
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D. Performing the Integrals 



Another way to compare the EFT results with the traditional ones is to carry out the 
integral Ii + I2 + I3 directly, and compare the final form of the resummed result. We wish 
also to show that the way we arrive at the final result is much simpler than the existing one 
in the literature. 

Specializing for the DY and Higgs case, the integral is then. 



h + h + h 



5 rf^ 



m 



We also need the solution of the renormalization group equation for as{fi^). Adopting the 
notation of Ref. lid we have 



1 



I bo 



In I 



+ 



I 



^(In^Z - In/ + / - 1) - -(In/ - r 
1% bo 



(81) 



where / = 1 + boasiQ"^) In fi'^/Q'^ and bi 
^(1) 



Ti+iPi- Let us start with the contribution of the 



(47r) 

term. Changing the integration variable from jj? to I, this contribution gives 




1 ^ri2N^iln/ 

1 -tts Q )-r-r 

bo I 



+ 



«s(Q' 



)'[| [ln^/-ln/ + /-l] 



4(ln/-l) 



21niV ^ ^ 



&oas(<5^ 



(82) 



where A = bQas{Q'^)\n.N . The last equation includes a pattern that repeats itself when 
other contributions are included. Taking as a working rule that InA^ ~ (l/as((5^)), the last 
two terms give rise to comparable contributions, however inside the curly brackets we have 
expansion in as{Q'^)- Thus the hierarchy is manifest. Carrying out the integrals in Eq. 
is very simple and we get 



'^1 



InA^ 



A 



(1) 
(9,9) 



47r6n 



2A + (1 -2A) ln(l -2A) 
A 



+ 



4> 



+as{Q^) 



2A + ln(l - 2A) + - ln2(l - 2A) 



A-nbi 



+ 2Aln(l - 2A) + - hi^il - 2A) 



1 - 2A 



(83) 



Expanding the A-terms in the last equation, we get a sum of the form ^^(Q^) ln"+^ from 
the first term, a"(Q^) In" A^ from the second term, and a""*"^ In" A^ from the last term. These 
are commonly called: leading logarithms (LL), next-to-leading logarithms (NLL) and next- 
to-next-to leading logarithms (NNLL), respectively. Higher logarithmic accuracies follow 
easily. 
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Consider now the contribution from A^^^y Similar to the A^^'^ contribution we get 



(47r)26c 



Jl-2X <■ 



X 21nA^ + 



/ - 1 



(84) 



(2) 

SO we see that A) s does not contribute to the LL but starts from NLL. This contribution 



IS 



A 



(2) 



^'^'^ [2A + ln(l - 2A)] - as{Q^)j^^ [2X + 2X^ + ln(l - 2A)] . (85) 



(47r)263 



(3) 

From the A, . term we get 



4(3) 

^(9,9) 



dl 



(4vr)36o A_2A P 
which is a NNLL contribution; 



-a^(Q2) [1 + Oiasm)] {2\nN + 



I - 1 

&0«s(<5^ 



(86) 



2\ {Q,a) 



2A 



(47r)362 1_2A' 



(87) 



The contribution from the term AB^^^ — Z*-*^ starts at NNLL accuracy since this term 
vanishes for i = 0,l. From Eq. §^ we have ABf^^^^- f^^^^ = {1/2){D'-^^ - ApoC2AfJ^g^). The 
contribution of this term gives 



1 



'S2 



(47r2)6oa.(g2^ 
which is a NNLL contribution; 

1 



'B2 



«.(g')- 



4/3oC2^ 



(1) 



D 



1-2A 



(2) 



{q,9) (<1,9) 



A 



(88) 



1-A 



(47r)26o 

Writing the sum of all contributions already obtained in the form of 

Ia, + Ia, + Ia, + Ib, = ^^Ng^^^^^ + ^^fj,) + «.(Q')^7fJ,), 



(89) 



(90) 
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we get 



A 



(1) 

(9,9) 



47r6n 



2A + (1 -2A) ln(l -2A) 



A 



A 



(47r 
4C2A 



(2) 

(^f^ [2A + ln(l-2A)] + 



^9,9)' 



(1) 

(9,S) 



2A + ln(l-2A) + -ln2(l-2A) 
2 



(2) 
(9,9) 



47r 
47r6|] 

^(9.g)^l 
(477)263 



(47r)26o 



1 - 2A 



2A + 2A ln(l - 2A) + - - 2A) 



2A + ln(l -2A) 



2A2 



1 - 2A 



2A 



(3) 
(9,9) 



(47r)362 i_2A 



[2A + 2A2 + ln(l-2A)] — 



2A' 



(91) 



The above functions sum the large logarithms to LL, NLL and N^LL, respectively. It is 
straightforward to get also the a'lg^'^^ which resumms the N^LL. It will contain contributions 



from ^l^^g) up to i 



4 and from D) . and D 



l(3) 



I (4) 



IS 



^^^^ uuxLyj. ^(^qgY The yet uncalculated quantity ^(^g) 
the only missing piece to complete the N^LL resummation program. The above results for 
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ii) 



agree with those in 0, i^]. We remind the reader that we have set the factorization 
scale and the renormalization scale equal to and the 7^; dependence is hidden in N used 
throughout. The analysis for the DIS case can be performed similarly and one also gets 
agreement with the known results. 



V. CONCLUSION 

Threshold resummation of logarithmic enhancements due to soft gluon radiation has been 
performed using the methodology of effective field theory. This method works to any desired 
(subleading) logarithmic accuracy and it is completely equivalent to the more conventional, 
factorization-based techniques. This has been illustrated to all three inclusive processes we 
considered: DIS, DY and the SM Higgs production. 

Conceptually and technically, however, this approach is much less complicated and it 
is physically more transparent than other ones. Working perturbatively in moment space 
(and for large values of N) we found that one does not need to introduce any additional 
nonperturbative quantities (other than the conventional PDFs), as is usually the case in the 
traditional approaches. All the quantities needed to get the resummed coefficient functions 
are straightforwardly obtained from fixed-order calculations of the form factors (which supply 
the C^*^ and the 7!*'*), the Altarelli-Parisi splitting kernels (which supply the 72*^) and the 
cross section for real gluon emission in the soft limit (from which we get the A4^'^). It should 
be mentioned that the given treatment of DIS is applicable only in the Bjorken limit where 
one takes to infinity first. However, for finite (but large) values of where the scale 
Q^(l — xY would emerge, a different treatment is needed. 

The method discussed in this paper can be extended straightforwardly to other processes. 
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their paper, the jet function is similar to the matching coefficient Ai here. 
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APPENDIX 



In this appendix, we collect the coefficient functions for deep-inelastic scattering, Drell- 
Yan and the Higgs production (within the large top-quark mass effective theory) to 0{a'^) 
in the soft limit of full QCD. They are used to extract the matching coefficients Ai in Eqs. 
(34-36). As we have remarked in the main paper, these results must be reproduced by 
calculations of an EFT in which only the soft and collinear degrees of freedom are taken 
into account. For DIS (see Refs. SHH) , Drell-Yan (see Refs. 0|) and Higg production 



see Refs. MM), we have 



16Pi(a;) + \TD^{x) + 5(1 - x) 



In^ 



+ 



241)2(2;) - VlV^{x) - (45 + 32(2)^^0(0;) 



+5(1 -x) 



51 



I2C2 + 400 



In 



-f8P3(x) - 18P2(a;) - (27 + 32(2)2^1(2;) + 2 x 48 x ( -- ) (3^^0(2;) 



+ (^ + 36(2 + 64(3 ) 1^0(2;) + 5(1 - a; 



331 



-5(1 - x) 
+5(1 -x) 
\CaCf 
+ 



^I?o(a;) + 5(l-x) 



InM ^ + 



+ 69(2 - 78(3 + 6(1 
8^ , , 58^ , , 



19 16 



-C2 



In 



4^ , ^ 58^ , , 
-V2{x) - -V^{x) 



247 8 



27 



C2 Vo{x) 



457 38 4 

+ + -(3 

36 3 3 



-f PoW- Y<5(l-i) 
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